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Abstract. The quenched Landau gauge gluon propagator is investigated in lattice QCD with large 
assimetric lattices, accessing momenta as low as q ~ 100 MeV or smaller. Our investigation focus 
/■"■n • on the IR limit of the gluon dressing function, testing the compatibility with recent solutions of the 
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Dyson-Schwinger equations. In particular, the low energy parameters k and a(0) are measured. 



^ . The confinement of quarks in hadrons and the chiral symmetry breaking mechanism 

are believed to be linked to the low energy properties of QCD. In particular, the gluon 
and ghost propagators can provide us information on the mechanism of confinement. 

Two first principles approaches to non-perturbative problems are Dyson-Schwinger 
equations (DSE) and lattice methods. In lattice QCD, the gluon propagator has been 
revisited a number of times. However, due to the lattice sizes used in previous studies, 
the access to the IR region was limited. For the Dyson-Schwinger solution for the 
gluon propagator, recently in the low energy behaviour of the propagator was 

solved analytically. Moreover, the authors found a parametrization for the gluon dressing 
function that fitted the numerical solution of the DSE. The lattice being complementary 
to DSE, allows a check of compatibility between the two methods. In particular, we 
would like to check for the behaviour of small and zero momenta gluon propagator. 

Our investigation uses pure gauge, Wilson action, SU(3) configurations, /3 = 6.0 
(a -1 = 1.885 GeV), on large assimetric lattices 1 : 16 3 x 128 and 16 3 x 256. They were 
generated using combinations of over-relaxation (OVR) and Cabibbo-Mariani (HB) 
updates. For the smaller (larger) lattice, a combined sweep of 7 OVR and 2 HB (7 
OVR and 4 HB) was used and 3000 combined sweeps for thermalization. The 160 (70) 
configurations were saved with a separation of 3000 (1500) combined sweeps. Note that, 
due to the large extension on the time direction, the lowest momenta considered is about 
93 MeV for the smallest lattice and 46 MeV for the largest lattice. 

For notation and details see [ 1 ] . The data for the Landau gauge gluon propagator 2 

<A»A>')> = V8(p + p>) S«*(S MV -^) ^p- (1) 
shows finite volume effects. To disentangle the finite volume effects, our data was 



1 All configurations were generated with MILC code |http : / /physics . indiana . edu/~sq/milc . html| 

2 For gauge fixing we used the steepest descent Fourier accelerated method described in 01 . 
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FIGURE 1. Gluon propagator. Statistical errors were computed using the jackknife procedure. 



compared with the results reported in [3]. No essential differences were seen between 
our pure temporal data and the Leinweber et al data. From now on, we refer only to 
pure temporal data - see figure [l] Note that the data seems to support a vanishing zero 
momenta gluon propagator D(p 2 ) . 

To test the compatibility between lattice and the DSE results [J], our data for the gluon 
dressing function was fitted to 
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-13/22, using two ansatz for the running coupling 
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with j8 = 11. 

The IR propagator was investigated fitting the lattice data to 
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Only when the first three lowest momenta of the larger lattice (\q\ < \39MeV) were 
considered, we were able to fit the IR analytical solution of DSE Zjri with an acceptable 
% 2 /d.o.f. = 0.39, meaning that the solution Z(q 2 ) ~ (q 2 ) 2K seems to be valid only 
for momenta below 150 MeV. The measured 3 k = 0.5003+29 ^ oes not P rov ide a clear 



Statistical errors for fit parameters were computed with 2000 bootstrap samples. 



answer about the zero momenta gluon propagator. The results for the IR2 approximation 
are shown in table UJ For both lattices, the fit to the largest range of momenta supports 
an infrared finite gluon propagator. 



TABLE 1. Zis2 fits. 



Lattice 




Range 


K 


A QCD (MeV) 


X 2 /d.o.f. 


16 3 x 128 
16 3 x 128 
16 3 x 256 


Ikl 
IM 
IM 


<461MeV 

< 553MeV 

< 644MeV 


0.5020+^ 
0.5122^2 

0.5199+.^ 


429+]° 
403+J 
377+5 


0.50 
1.41 
1.19 



The fits to all data distinguishes the functional forms CCi(p 2 ) and 0C2(p 2 ) for the 
running coupling constant. For <X\(p 2 ), the % 2 jd.o.f. > 2. Moreover, a(0) ~ 10 is 
much larger than the DSE result: cc(0) = 2.972 JHIH]. However, if one uses C^Cp 2 ), 
the lattice data is well described by © - see table |2j oc(0) can be computed from the 
asymptotic behaviour of the QCD /3 function and 0C2(p 2 )- Then, a(0) = (47r//3o)a2 
giving a(0) = 2.78^2 an d 2.74+2 f° r the smaller and larger lattices, respectively. The 
analysis of all momenta favours a finite zero momenta gluon propagator. Our figures for 
K are smaller than those reported in the DSE studies, K ~ 0.595 04J]. However, our K 
values reported in tableEJare within the figures discussed in |0,[li. 



TABLE 2. Fits to all lattice data. 



Lattice 


K 


Aq CD (MeV) 


a\ 




a 2 


X 2 /d.o.f. 


16 3 x 128 
16 3 x 256 


0.5439+ 3 ^ 
0.5314+25 


352+J 
354^ 3 


0.0063; 
0.0065; 


f4 
-3 
f4 

-3 


2.43+ 2 
2.40^ 


1.74 
1.56 



We are currently involved in improving the statistics of our analysis and, simultane- 
ously, trying to understand the role of Gribov copies 4 relying on the method discussed 
in [9]. We plan also to compute the ghost propagator and running coupling constant 
directly from the lattice. 
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4 In CI > tn e effect of the Gribov copies was estimated as a 2 to 3 (J effect. This does not change our 
predictions for the zero momenta gluon propagator. 



